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Normal forms for pseudo-Riemannian 2-dimensional metrics whose 
geodesic flows admit integrals quadratic in momenta 

Alexey V. Bolsinov* Vladimir S. Matveevj Giuseppe Pucacco* 



Abstract 

We discuss pseudo-Riemannian metrics on 2-dimensional manifolds such that the geodesic flow admits 
a nontrivial integral quadratic in velocities. We construct (Theorem [TJ local normal forms of such 
metrics. We show that these metrics have certain useful properties similar to those of Riemannian 
Liouville metrics, namely: 

• they admit geodesically equivalent metrics (Theorem [2}; 

• one can use them to construct a large family of natural systems admitting integrals quadratic in 
momenta (Theorem HJ; 

• the integrability of such systems can be generalized to the quantum setting (Theorem [5}; 

• these natural systems are integrable by quadratures (Section 12. 2. 5) . 
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1 Introduction 

Consider a pseudo-Riemannian metric g = (gij) on a surface M 2 . A function F : T*M — > M is called an 
integral of the geodesic flow of g, if {H, F} = 0, where H := \g l ->PiPj : T*M — > R is the kinetic energy 
corresponding to the metric. Geometrically, this condition means that the function is constant on the orbits 
of the Hamiltonian system with the Hamiltonian H . We say the integral F is quadratic in momenta if, 
in every local coordinate system (x,y) on M 2 , it has the form 



a(x, y)p\ + b(x, y)p x p v + c(x, y)p 2 , 



(1) 



with (x,y,p x ,p y ) canonical coordinates on T*M 2 . Geometrically, formula |1| means that the restriction of 
the integral to every cotangent space T*M 2 = M 2 is a homogeneous quadratic function. Of course, H itself 
is an integral quadratic in momenta for g. We will say that the integral F is nontrivial, if F ^ const • H 
for all const € M. 

The main result of this paper is Theorem [T] below, which gives us a list of local normal forms of metrics of 
signature (+, — ) whose geodesic flows admit a nontrivial integral quadratic in momenta. For the Riemannian 
case (and, therefore, for the signature (— , — )) such metrics are the well-known Liouville metrics. 

Theorem 1. Suppose the metric g of signature (+, — ) on M 2 admits a nontrivial integral quadratic in 
momenta. Then, in a neighbourhood of almost every point there exist coordinates x, y such that the metric 
and the integral are as in the following table: 





Liouville case 


Complex- Liouville case 


Jordan-block case 


9 


(X(x)-Y(y))(dx 2 -dy 2 ) 


^s{h)dxdy 


(l + xY'(y))dxdy 


F 


X(x)-Y(y) 


2 2 , n^( h ) 
Px-Py + 2 WUT)PxPy 


Px 2 l+xY'(y)P x Py 



where 5R(/i) and 3(/i) are the real and imaginary parts of a holomorphic function h of the variable z := x+i-y. 
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Given a metric and the quadratic integral, it is easy to understand what case they belong to. Indeed, 
for the integral (|TJ) the matrix 

Fi, =(l !) 

can be viewed as a (2, 0)-tensor: if we change the coordinate system and rewrite the function F in the new 
coordinates, the matrix changes according to the tensor rule. Then, 

Gj (2) 

a 

is a (1, l)-tensor. By direct calculation we see that Gj has two different real eigenvalues in the first case, 
two complex-conjugate eigenvalues in the second case and is (conjugate to) a Jordan-block in the third case. 
This also explains our choice of the names for the normal forms of the metrics. Indeed, in the Riemannian 
case, the tensor @ always has two real eigenvalues. In particular, the normal form of the Riemannian 
metric admitting an integral quadratic in momenta, which is traditionally called Liouville form (or Liouville 
metric), is very similar to the metric of our "Liouville" case. One can view our "Complex-Liouville" case 
as the complexification of the standard Liouville metric: if in the expression 

(X(x)-Y(y))(dx 2 +dy 2 ) 

we replace X by (a holomorphic function) h{z), Y by h(z), dx by dz, and dy by idz, we obtain the Complex- 
Liouville metric up to the factor 8i. The Jordan-block case has no direct analog in the Riemannian setting. 

Remark 1. The corresponding natural Hamiltonian problem on the hyperbolic plane has recently been 
treated in [39] following an approach used by Rosquist and Uggla [40] . Systems with indefinite signature 
have been investigated before in the classical works by Kalnins and Miller on separation of variables [20l [38] , 
see also [T2] [TCJ [35] . Other possible approaches are based on Killing tensor theory [5] , r- matrix theory [53] 
and algebraic methods [15] . For the corresponding quantum case we refer to |19j and references therein. 

Remark 2. A part, if not all credits for the results of the present paper should be given to Darboux, see 
[14], §§592-594,600-608]. There is no doubt that Darboux was very close to Theorem [TJ to the results of 
Section 12.2.21 and, to a certain extent, to Theorem [5] of our paper, and could get it if he would have been 
interested in the pseudo-Riemannain metrics. More precisely, 

• In [141 §593], Darboux gets the Riemannian Liouville metrics. Since he worked over complex coordi- 
nates, his formulas can be interpreted as our Liouville and Complex-Liouville cases. 

• In [14, §594], Darboux gets (a case that could be interpreted as) the Jordan-block case. 

• The formulas of Section [2.2.21 of the present paper are similar to that of [131 §594]. 

However, Darboux was interested in the positive definite metrics only. Actually, in his time it was unusual 
to consider indefinite metrics, since the applications of pseudo- Riemannian metrics to general relativity and 
cosmology appeared much later. Darboux worked over complex coordinates x, y and explicitly remarked on 
the transformation x — u + iv,y — u — iv leading to the standard metric of the (+,+) case, with no mention 
of a possible interpretation of x, y as real coordinates. The only exception is the Jordan- block case with 
constant function Y (equations (24,25) of [HI §594]), where one can get the surfaces of revolution. 

The results of this paper were announced in [10] . 

2 Applications 

2.1 Applications in geometry: normal forms for 2-dimensional geodesically 
equivalent metrics 

Two metrics g and g on one manifold arc geodesically equivalent, if every (unparametrized) geodesic 
of the first metric is a geodesic of the second metric. Investigation of geodesically equivalent metrics is a 
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classical problem in differential geometry, see the surveys [3, 33, 37J or/and the introductions to [3"Tl |3"21 151] . 
In particular, normal forms for geodesically equivalent Riemannian 2-dimensional metrics were already 
constructed by Dini [17j . An easy corollary of Theorem [T] is the following theorem which gives normal forms 
of geodesically equivalent nonproportional metrics such that one of them has signature (+, — ). 

Theorem 2. Let g, g be geodesically equivalent metrics on M 2 such that g has signature (+,—), and 
g =/= const ■ g for every const G K. Then, in a neighbourhood of almost every point, there exist coordinates 
such that metrics are as in the following table: 



Liouville case 



Complex-Liouville case 



Jordan-block case 



(X(x)-Y(y))(dx 2 -dy 2 ) 



^s(h)dxdy 



(1 + xY'(y))dxdy 



( 1 1 A ( dx 2 _ dy 2 

{YJyj Xffi) [x{Fj Y(y) 



( 9W V rf _a 
\S(h) 2 +Si(h) 2 J "■' 



2t^ 



Sft(fe)S(ft) 



r 3(h) 2 +}R(h) 2 



dxdy 
''dy 2 



l^l(-2Y(y)dxdy 
+ (l + xY'(y))dy 2 ) 



where h is holomorphic function of the variable z := x + i ■ y . 

Remark 3. It it natural to consider the metrics from the Complex-Liouville case as the complexification of 
the metrics from the Liouville case: indeed, in the complex coordinates z = x + i- y,z = x~i-y, the metrics 
have the form 



ds] = -\{h(z) - h{z)) (dz 2 - dz 2 ) 

,7„2 _ 1 \ ( dz 2 dz 2 ' 

b 9 4 h(z) J h(z) 

Remark 4. In the Jordan-block case, if dY ^ (which is always the case at almost every point, if the 
restriction of g to any neighborhood does not admit a Killing vector field), after a local coordinate change, 
the metrics g and g have the form (see also Remark [T4"|) 



ds 2 g = 



ds 2 = 



Y{y) 
2(Y(y) 



xj dxdy 



-dxdy 



{Y(y) + xf 
y 4 



dy 2 



Proof of Theorem [U We will use the next theorem which probably was already known to Darboux [14j 
§608] . For recent proofs, see [Ml [571 Hi H • 

Theorem 3. Let g be a metric on M 2 and h € F(S f 2M 2 ) be a symmetric nondegenerate bilinear form on 
M 2 . Consider the following metric 



/det( fl )V 



(3) 



y \det{h) J 

on M 2 . If g and g are geodesically equivalent, then the function 

h : TM — > R, h(£) :=h(£,£) 
is an integral for the geodesic flow of g. 

Remark 5. Theorem [3] and Corollary [T] below bear some resemblance with other classes of transformations 
between dynamical systems [TJ [21 [TTJ [38J STJ [42J [43]. However, the present result is of different nature and 
is deeper because, in order to construct the second system, one needs to know the quadratic integral of the 
first one. 

Combining Theorem [3] with Theorem [T] we obtain that, in a neighbourhood of almost every point, 
geodesically equivalent metrics g and g are as in the table in Theorem [2] (we assume that g has signature 
(+, — ) and that g ^ const-g). Thus, in order to prove TheoremfSl we need to show that the metrics from the 
table are indeed geodesically equivalent, which can be done by direct calculations. Indeed, it is well-known, 
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see for example [T5| §40 of Ch. Ill], that two metrics are geodesically equivalent if and only if the difference 
of their Levi-Civita connections has the form Tj5\ + Yfc£>] for a one-form T = (T^). Direct calculation of 
the Levi-Civita connections for the metrics shows that it is indeed the case: the form T equals 



2 V X(x) Y(y) 
for the normal forms of the metrics in the Liouville case, 

%(h)JL%(h) + Vi(h)^(h) ^ ^(h)^(h)-^t(h)^(h) d 

(%(h)) 2 + mh)) 2 x+ (z(h)f + (M(h)f v 

for the complex Liouville case and yr^ dy for the Jordan-block case. □ 

Corollary 1. Let g be a metric on M 2 and h S r(«S , 2-M 2 ) be a symmetric nondegenerate bilinear form on 
M 2 . Then, g and the metric ([3]) are geodesically equivalent, if and only if the function 

k : TM -> K, h(£) = 

is an integral for the geodesic flow of g. 

Proof. In the direction the statement coincides with Theorem [3J In order to prove in "<^=" 

direction, it is sufficient to check the statement in the neighbourhood of almost every point. Here, the 
metrics g, g and the integrals h are given by Theorems 11121 and are related precisely by formula ([3]). □ 

Remark 6. Theorem [3] had found a recent important application in the solution of two problems explicitly 
stated by Sophus Lie in [55] due to [T2] [35] . 



2.2 Applications in mathematical physics 

2.2.1 Natural systems admitting an integral quadratic in momenta 

For a pseudo-Riemannian manifold (M,g), a natural Hamiltonian system is a Hamitonian system with 
H : T* M — > R of the form H := H g + U = hg^PiPj + U(x,y). We say that a natural Hamiltonian system 
is quadratically integrable, if there exists a function F of the form F = F g + V = F^piPj + V(x, y) such 
that {H, F} = with F ^ consti • H + const2 for all consti, const2 £ M. 

Remark 7. In [39], the natural Hamiltonian system on the Minkowski plane has been reduced to the corre- 
sponding kinetic Hamiltonian system with conformal (Jacobi) pseudo-Euclidean metric. 

Theorem 4. Let g be a metric of signature (+, — ) on M 2 . Assume a natural Hamiltonian system with 
Hamiltonian H g + U to be quadratically integrable with integral F = F g + V. Then, in a neighbourhood of 
almost every point, there exists a coordinate system such that the metric g and the functions F g , U , V are 
as in the following table: 





Liouville case 


Complex- Liouville case 


Jordan-block case 


9 


(X(x)-Y(y))(dx 2 -dy 2 ) 


^{h)dxdy 


(1 + xY'(y))dxdy 


F g 


X(x)pi ) -Y(y)pi 
X(x)-Y(y) 


Px Py ' *<5(h)P x Py 


Px A \+xY'(y)P x Py 


U 


1 X{x)-Y(y) 
2X(x)-Y{y) 


z(h) 


xY;(y)+Y 2 (y) 
l+xY'(y) 


V 


Y(y)X(x)-X(x)Y(y) 




Y xY;(y)+Y 2 (y) , y l \ 


X(x)-Y(y) 




1 l+xY'(y) +I H!/J 



where h, hi are holomorphic functions of the variable z := x + i ■ y . 

Proof. It is well known (see, for example, [5]), that the condition {H, F} = is in this case equivalent 
to the following two conditions: 

{H B ,F g } - 0, (4) 
2dUoG = dV, (5) 
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where G is given by @. In tensor index notations, ^ is 



dU dV 

2G*™ = ™. (6) 
J dx 1 dxi K ' 

Indeed, condition {H, F} = is equivalent to the following equation: 

{H g ,F g } + {H g ,V}-{F g ,U} = 0. 

Since {H g , F g } (respectively, {H g , V} — {F g , U}) is a third degree-polynomial in momenta (respectively, first 
degree), the latter equation is equivalent to: 

{H g ,F g } = (7) 
{F g ,U} = {H g ,V}. (8) 

We see that Q coincides with (Q} and ([5]) is equivalent to 

nT TiidU udV 
2F 3 r = g 3 r , 

dx l dx % ' 

which is equivalent to ([6]) and therefore to (JSJ). 

Condition Q tells us that the function F g is an integral quadratic in momenta for the geodesic flow of 
g. Clearly, F g is nontrivial. Indeed, if F g — consti • H g , then condition (JSJ) reads consti o dU = dV implying 
V = consti • U + const2. These in turn imply F = consti • H + const2, which contradicts the assumptions. 

Thus, F g is a nontrivial integral of the geodesic flow of the metric g. By Theorem [TJ almost every point 
has a neighbourhood with local coordinates (x, y) such that g and F g are as in the table. In order to prove 
Theorem |4l it is sufficient to show that, for every column of the table, the functions U and V are complete 
solutions of equation ([5]). Here we consider the three cases in detail. 

Liouville case. Assume g, F g are as in the first column of the table. Then the form dU o G is 

-y { y)f x dx-X { x) d ^dy 



and condition (HI) reads 



f dY(y)U _1§V 

s ox(x)u _i§v 

{ dy 2 dy ' 

Differentiating the second equation w.r.t. x and subtracting the derivative of the first equation w.r.t. y, we 
obtain 

n 9 (tt! \ dU \ 9 fvt ^ 9U \ d*(X(x)-Y(y))U 
° = to [ X{X) dy- ) { Y{y) to J = dxTy 

implying 

l X{x)-Y{y) 
2 X(x)-Y(y) 

for certain functions X = X(x) and Y = Y(y). Substituting U in ([9]), we obtain 

X(x)Y(y) - Y(y)X(x) 



V = 



X(x)-Y(y) 



Thus, in the Liouville case, U and V are as in the table. 

Complex-Liouville case. In this case 2dU o G is equal to 



f d^{h)U _ dSt(h)U \ ^ f d^{h)U &X(h)U \ 
\ dx dy J \ dy dx J 



dy 
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and condition is equivalent to the following system of PDE: 

{ ddt(h)U _ dQ(h)U _ dV 
dx dy dx ' ( AC\\ 

d?R(h)U d<3(h)U _ OV { W > 
dy dx dy ' 

We see that these equation are precisely the Cauchy-Riemann condition for the function hi := $t(h)U — V + 
i ■ %(h)U. Thus, 

rj _ gfo) 
9f(ft) 

and 

V = mh)U - Mhi) = KW^rr - ^(hi) ■ 

We see that U and V are as in the table. 

Jordan-block case. In this case the 1-form 2dU o G is 

-Y(y)^dx + ((1 + xY'(y))^ - Y{y)^\ dy 
and condition ((SJ is equivalent to the following system of PDE: 



(11) 



The first equation in (fTTj) is equivalent to V = —Y(y)U + Y\(y). Substituting this in the second equation, 
we obtain 

,dU dY(y)U . 
(l + X Y> { y))--Y { y)- = -^ + Y{ { y) 

which implies 

ajujpm _ Y , {v) 

and therefore (1 + xY'(y))U = xY{(y) + Y 2 (y). Thus, 

xY{(y)+Y 2 (y) 
l + xY'(y) 



□ 



2.2.2 Integration by quadratures of natural systems admitting an integral quadratic in 
momenta 

Since the time of Jacobi it is known that (in the 2-dimcnsional Ricmannian case) nontrivial integrals 
quadratic in momenta are extremely helpful for the description of dynamics of natural systems: indeed, 
in this case 

• the Hamilton equations, which are a system of four ODE on T*M 2 , can be reduced to a parameter- 
depending system of two ODE on M 2 . 

• Moreover, it is possible to construct a characteristic (= function constant on the solutions) of this 
system by means of the integration of certain functions of one variable only. 



G 



See E5] for details. 

Classically, the second property is referred to as "the system is integrable by quadratures". Both 
properties are useful for exact solutions, for numerical analysis and for a qualitative description of (the 
solutions of) the Hamilton equations. We are going to show that these nice properties persist in the pseudo- 
Riemannian setting. 

Liouville case. There is virtually no difference with respect to the Riemannian setting. Consider 
H = H g + U and F = F g + V such that g, F g ,U,V are as in the first column of the table from Theorem [H 
Then, the first two Hamilton equations are 



_ dIL — 
dtV = dp y ~ X-Y 



it dp* X-Y ' f-i \ 

d _ an _ p y U^j 



Since the functions F and H are constant on the solutions of the system, for every point (x, y,p X iPy) of 
the solution we have 



1 X(x)-Y(y) 



2 X(x)-Y(y) ^ 2 X{x)-Y{y) 
X(x)p 2 y -Y{y)pl Y{y)X(x)-X(x)Y{y) _ 



Ho j 

Fn. 



X(x)-Y(y) + X(x)-Y(y) r ° 

This is a linear system on p%,Py, solving it w.r.t. p x and p y we obtain 

pi = 2H X(x)+F -X(x), 
V\ = 2H Y(y)+F -Y(y). 



(13) 



Substituting these in (fT2|) . we obtain 

d „ 



dt" ~ £l x-y '■— Vl ' Q4\ 

d = V2H Y{y)+F -Y{y 1 

dt y t 2 X _ Y ■ vi. 

We see that Hamilton equations can be reduced to a system of two ODE on M 2 depending on the parameters 
#0,^0 S K and e t G {-1,+1}. 

Clearly, a function K(x,y) is a characteristic of the system (fT4|) if dK vanishes on the vector field 
v := (vi,V2). Since the form 

sidx e 2 dy 



'2H X(x) + F - X(x) \f2H Y(y) + F - Y(y) 
vanishes on v and is closed, the function 

m :=( r B=r ^ - ^ 



I'll 



2H X(Z) + F - X(0 J vo J2H Q Y(0 + F - Y(£) 



is a characteristic. We see that in order to find a characteristic, we only need to integrate two functions of 
one variable each, i.e., the system is integrable by quadratures. 

Complex-Liouville case. Consider H = H g + U and F = F g + V such that g, F g ,U,V are as in the 
second column of the table from Theorem dJ Then, the first two Hamilton equations are 

J dt X dp* 3(h)' /ir^ 

1 ±n, - M. - 2 p* ^ > 

y dty ~ dp y ~ Wjij ■ 

Since the functions F and H are constant on the solutions of the system, for every point (x,y,p x ,p y ) of 
the solution we have 



Q PxPy | Will _ jl 



ri-^ + R W (2^ + ^)-R(fc 1 ) = F . 
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Subtracting the first equation times from the second, we obtain 

vl-vl = -(»(ft)fr -K(fci)) + Jb- 

From these, adding (respectively, substracting) to (respectively, from) the second equation the first equation 
times i, we obtain 

{Px-i-Pyf = -{H 3t(h)-to{hi)-F )-i-(H,&(h)-Q{hi)) = -Hoh + hi + Fo, 
(p x + i-Pyf = -W^-W-^+i-W)-*)) = -H h + hi+F . 

Remark 8. Since \{p x — i ■ Py) is the canonical momentum conjugate to z = x + i ■ y, these equations are 
the complex analog of (|13p. 

Then, = eSft (y/^Hoh + h\ + F ) and p y = — e3 (^/—H^h + h[ + F ) (the choice of the branch of the 
square root is hidden in e). Substituting these in (|15[) . we obtain 

_d _ -2ga(v/-gofe+/ti+Fo) _ 

dt x ~ ~ : ~ "ii / lfi s 

d 2eK(V--H"o?i+/ii+F ) v ' 

= — — *m : = V2 - 

We see that Hamilton equations can be reduced to a system of two ODE on M 2 depending on the parameters 
H Q ,F Q G E, and £ e {-1,+1}. 
Consider the 1-form 

9ft(V-ffofe + /n+F ) g (V-ffpfe + fei + g) 
' \-H h + hi + F \ + \-H h + hi+F \ V - 



The Cauchy-Riemann conditions for the holomorphic function \J—H$h + hi + Fq imply that the form is 
closed. Clearly, the form vanishes on the vector field v = (v\,V2)- Then, the function 



Hoh + hr + Fol J yo \-H h + hi+F { 



is constant on the solutions of (|16p. i.e., is a characteristic of the system. It is easy to check by direct 
calculations that in the complex coordinate z the form B is 



25R 

Thus, the function K equals to 

25R 



dz 



\/—Hoh + h 1 +F Q 



^-H h(0 + h 1 (0+F 0/ 



i.e., the system is integrable by quadratures. 

Jordan-block case. Consider H = H g + U and F — F g + V such that g, F g , U, V are as in the third 
column of the table from Theorem 0J Then, the first two Hamilton equations are 

( At = ML = 2p y 

J dt X dp* l+xY'(y)> Mn\ 

1 A, , = ML 2 Px K 11 ) 

{ dtU dp v l+xY'(y) ■ 

Since the functions F and H are constant on the solutions of the system, for every point (x, y,p X iPy) of 
the solution we have 

o PxPb i Y 2 (y)+xY((y) _ n 

^l+xY'iy)^ l+xY'(y) ~ 170 » 



s 



Adding the first equation times Y(y) to the second one, we obtain 



pi = H Y(y)-Y 1 (y) + F Q 

2 PxPy = x(H Y'(y)-Y{(y))+H -Y 2 (y) 



where e E {—!,+!}. Substituting these in (jT7j), we obtain 




e y /H Y(y)-Y 1 (y)+F , 

, x(H Y'(y)-Y{(y)) + H -Y 2 (,y) 
y / H Y(y)-Y 1 (y)+F 



x(ff Y'(v)-Y({y))+H -Y 2 { V ) 

—x = e , ; / , == — = := vi , 



y 



' (i+gr' fa))VH yfa)-yi(i/ )+fb 

2^H Q Y{y)-Y 1 {y)+Fo 
■ l+xY'(y) 



(18) 



v% ■ 



We see that Hamilton equations can be reduced to a system of two ODE on M 2 depending on the parameters 
H Q ,F G E, and e E {-1,+1}. 
Consider the 1-form 



D 



dx 



lx(H Y / (y)-Y{(y))-Y 2 (y) + H 



^H Y(y)-Y 1 (y) + F 



^H Q Y(y)-Y 1 (y)+F 



(H Q Y(y)-Y 1 (y) + F f 2 
1 Y 2 (y)-H 
^(H Q Y(y)-Y 1 (y) + F f 2 



dy 



dy. 



(19) 
(20) 



By (|2"01) , the form is closed. By (fH)|) . the form vanishes on the vector field v — {vi,v 2 ). Then, the 
function 



K(p) 



D 



^/F Q -Y 1 (y) + H Q Y(y) 



Y 2 (Z)-H 



(Fo-Y^+HoYiOf 2 



d£ 



is a characteristic of the system (jTSJ) , i.e. the system is integrable by quadratures. 
2.2.3 Quantum integrability 

Let g be a metric, and (F v ) E T(S 2 M 2 ) be a symmetric bilinear 2-form on T*M 2 . Consider the following 
two linear partial differential operators A g ,F g : C°° — + C°°: 



^ V|de*(ff)l 9x * dx J 

Remark 9. The first operator is the Beltrami-Laplace operator of the metric g; another way to write it down 
is 

i,3 

where V is the Levi-Civita connection of g. The second operator is a natural quantization of the function 
X)i j F^PiPj and another way to write it down is 



In particular, both operators do not depend on the choice of the coordinate system. 

Remark 10. The symbols of A g and of T g are — 2H := —2 Yli j Q^PiPj an d F^PiPj: respectively 
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Theorem 5. Let F = ^\ . F lJ 'piPj + V(x,y) be a quadratic integral of the natural Hamiltonian syste 
5 Si j d^PiPj + U(x,y) on T* M 2 . Then, the operators 

H := Ag - 2U 

and 



T g + V 



commute: TL o T = T o H. 



Remark 11. The Riemannian analog of Theorem [5] follows from 9 , 23, 2~9"1 130j. 

Proof of Theorem [H It is sufficient to check the statement at almost every point, i.e., for the metrics 
and the integrals from Theorem [U Direct calculations shows that in this case the operators A g and T g are 
as in the following table: 





Liouville case 


Complex- Liouville case 


Jordan-block case 






-1 ( d ' 2 d 2 \ 


-4 d 2 


-4 d 2 


X(x) 


-Y(y) \dx^ dy 1 ) 


^s(h) dxdy 


l+x 1 Y'(y) dxdy 




1 


(X(x)^-Y(y)£) 


d 2 d 2 | o «C0 d 2 


d 2 o Y(y) 2 


X(x)-Y(y) 


dx 2 dy 2 ^(^) dxdy 


dx 2 l-\-xY' (y) dxdy 



where h is a holomorphic function of z = x + i ■ y. 



To prove that H = A g 
Jordan-block cases: 



2U and T = T a + V commute, we first observe that in the Liouville and 



T g + V 



dx 2 



f ■ (A g - 2U) + h, 



where / = X(x), f\ — X(x) for the Liouville case, and / = and fx = Yi{y) for the Jordan block case. 
Similarly, in the complex Liouville case, we have 



Fa 



d 2 d 2 



where / = -^,/ 1 = «(M 

The Laplace-Bcltrami operator A g in all the cases is of the form A g = A 1 A go , where A ff0 is the 
Laplace-Beltrami operator of the flat metric go (more specifically, go is dx 2 — dy 2 in the Liouville case, and 
2dxdy in the complex Liouville and Jordan-block cases). Using the fact that A go commutes with J^, it is 
straightforward to verify the following commutator formula: 



[A g - 2U, 



— 

dx 2 



^xx 



\ x d 
> t 

A dx 



{A g -2U) + 2 



(XU)a 

A 



+ 4 



(XU)x d 
A dx 



(here we use standard notation for the commutator of two linear operators [A, B] = Ao B — B o A). 
The two following formulas are standard: 



[A g - 2U, f ■ (A B - 2U)] = (A g f - 2grad fl /) o (A g - 2U) 



and 



[A ff -2C/,/ 1 ]=A s / 1 -2grad ff / 1 , 

where the vector field grad g / is viewed as a first order differential operator, i.e., grad ff / = <7*" ? J^r gfj ■ 
Thus, in the Liouville and Jordan-block cases, we have: 



[H,7] = [A g -2U,T g + V] 



(XU)a 

A 



2^J-+A 9 /-2grad 9 /) o(A 9 -2('i 



A dx 
[ {XU)xd_ 
A dx 



AJ 1 -2grad / 1 
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Hence, the commutativity condition \H 1 T] = 7i o T — 7Y o J 7 = splits into four simple equations (here we 
use the fact that A g = A _1 A fl0 and grad g = A _1 grad fl0 ): 

« ^^-grad 9o / = 

(ii) A xx + A 9o / = 

(hi) 2(AC/)^-grad go / 1 = 

(iv) 2(AE/) xx + A 5o / 1 =0 

Each of these equations has natural meaning. Indeed, (i) and (ii) mean that the operators A g and T g 
commute (without potentials), (iii) and (iv) give the "new" commutativity conditions involving the poten- 
tials. The first and third equations are equivalent to the commutativity of classical integrals, whereas the 
second and the fourth keep additional "quantum" information. It is interesting to notice that the quantum 
conditions (ii) and (iv) can be obtained from the classical ones (i) and (iii) by "differentiating" so that in 
our particular case the quantum integrability in dimension 2 turns out to be a corollary of the classical one: 

\ xx + A go f = div(A x ^ - grad 90 /) 

2(\U) XX + A go h = div(2(AC/)^ - grad go /x) 

However, each of the above four conditions can be verified directly. Taking into account the following 
explicit formulas: 

Liouville case : 

A 9o =-^ + ^, S rad ff o/ = /^-/^. X = X(x)-Y(y) 
Jordan-block case : 

A 9 o=-2af^ l?^*/ = + * = 1/2(1 + xY'(y)) 

we see that equations(i)-(iv) become: 

Liouville case: Jordan-block case: 

(A x — fx)g^ + fy-Qy- = (A x — fyjfa —fxQy- =0 

^xx fxx ~t~ fyy X xx 2f X y 

(2(XU) X - (A) x )£ + {h) y -§- y = (2(XU) X - {h) y )4L - ( h ) x -§- = 

2(XU) XX - {h) xx + (fl)yy = 2{XU) XX - 2{h) xy = 

and obviously hold for A, / and fx indicated above. 

The complex Liouville case is absolutely similar, the only difference is the additional term -^p, which 
leads to the following system of relations: 



Complex Liouville case: 

f >~5v~ 



(A x - fy)-§£ + (-A y - /a;)gjj — 
X-xx ^yy 2f X y 

(2(XU) X - {fx) v )-§- x + {-2{XU) y - (h) x )-§- y = 



2(A!7) XX - 2{XU) yy - 2{h) xy = 
each of which obviously holds for / = fx = A = 2XU = Qf(/i 2 ). □ 

3 Proof of Theorem [I] 

3.1 Admissible coordinate systems and Birkhoff-Kolokoltsov forms 

Let g be a pseudo-Riemannian metric on M 2 of signature (+, — ). Consider (and fix) two vector fields Vx, V% 
on M 2 such that 

• g(Vx,Vx) = g(V 2 ,V 2 ) = and 
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• g{V u V 2 ) > 0. 

Such vector fields always exist locally, (and since our result is local, this is sufficient for our proof). For 
possible further use, let us note that such vector fields always exist on a finite (at most, 4-sheet-) cover of 
M 2 . 

We will say that a local coordinate system (x, y) is admissible, if the vector fields and ^ are 
proportional to Vi, V 2 with positive coefficient of proportionality: 

d d 

— =\ 1 {x,y)V 1 (x,y), — = X 2 (x,y)V 2 (x,y), where A* > 0. 

Obviously, 

• admissible coordinates exist in a sufficiently small neighbourhood of every point, 

• the metric g in admissible coordinates has the form 

ds 2 = f(x, y)dxdy, where / > 0, (21) 

• two admissible coordinate systems in one neighbourhood are connected by 

%new 1 I X ne w{Xold)\ , dx new dy new fnn\ 

= ) ( i where — > 0, — > 0. (22) 

VnewJ \ynew[yold) J dx ld dy a ld 

Lemma 1. Let (x,y) be an admissible coordinate system for g. Let F given by ([T]) be an integral for g. 
Then, 

Bl := h ] ^ dx > I respectively, B 2 := ] dy ] 

is a 1-form, which is defined at points such that n ^ (respectively, c ^ 0). Moreover, the coefficient a 
(respectively, c) depends only on x (respectively, y), which in particular implies that the forms B\, B 2 are 
closed. 

Remark 12. The forms B\,B 2 are not the direct analog of the "Birkhoff" 2-form introduced by Kolokoltsov 
in [22]. In a certain sense, they are a real analog of the square root of the Birkhoff form. 

Proof of Lemma [TJ The first part of the statement, namely that 

zdx, [ respectively, — dy | 



\/Hx,y)\ \ ' ' yj\c(x,y)\ 

transforms as a 1-form under admissible coordinate changes is evident: indeed, after the coordinate change 
([22j). the momenta transform as follows: p Xold = Px nau , d dx"Z 1 P x °id ~ Pa„ em ■ Then, the integral F in 
the new coordinates has the form 

dx new \ 2 . dx new dy new ( dy new \ o 

a Px new + ~j ~j Op Xnew Py new + — CP 

dxoid J dxoid dy id \ d Vold J 



Then, the formal expression } dx id (respectively, } dy id) transforms into 

v/M V l<=| 



1 dx oM / 1 dy oM 

dx new respectively, — =- dy ne 



\a\ dx new y ^ \/\c\dy 

which is precisely the transformation law of 1-forms 
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Let us prove that the coefficient a (respectively, c) depends only on x (respectively, y), which in particular 
implies that the forms B\, B 2 are closed. If g is given by (j2"T]) , its Hamiltonian is 

ZPxPy 



H 



and the condition {H, F} — reads 



J 2p x Py 2,1, 1 2 



2 

i.e., is equivalent to the following system of PDE: 



j2 \PUf a v) + P 2 xPy(fa x + fb y + 2f x a + f y b) + p y p x (fb x + fc y + f x b + 2/„c) + P y (c x f)) 



fa x + fby + 2f x a + f y b = , 



a y = , 

(23) 



fb X +fCy+f X b+2fyC = 0, 

c x = 0. 

Thus, a = a(x), c = c(y), which is equivalent to state that B\ :— ) dx and B2 ■= —j=dy are closed 

VH vl c l 



forms (assuming a ^ and c 7^ 0). □ 
Remark 13. For further use let us formulate one more consequence of equations (|23[) : if a = c = in a 
neighbourhood of a point, then bf = const, implying F = const • H in the neighbourhood. 

Assume a 7^ (respectively, c 7^ 0) at a point pq. For every pi in a small neighbourhood i7 of po consider 

/ \ 

(24) 



7 : [0, 1] U 



respectively, y new := J B 2 

V 7 : [0, 1] -> C/ 



with 7(0) =po,7(l) =pi. 

Locally, in the admissible coordinates, the functions x neu , and j/net« are given by 

f x 1 /" y 1 

a;neiu(x) = / 7= rff, y n ew(y) = / , , dt. (25) 



The coordinates (x new ,y ld), ((% old, ynew), (x new ,y new ), respectively) are admissible. In these coordi- 
nates the forms Bi,B 2 are given by dx new , dy new implying that a = c = ±1 (more precisely: a new = 
sign(a o;d ), c new = sign(c oW )). 

3.2 Proof of Theorem ffl 

We assume that g on M 2 of signature (+,-) admits a nontrivial quadratic integral F given by (Q}. Consider 
the (1, l)-tensor G given by ([2]). In a neighbourhood of almost every point, the Jordan normal form of this 
(1, l)-tensor is one of the following: 

Case 1 f q ^) j wne re A, /i 6 E. 

Case 2 ^ ^ ^ , where A, /i 6 E. 

Case 3 j^V where A £ E. 

Moreover, in view of Remark I13[ there exists a neighbourhood of almost every point such that A ^ /1 
in case 1 and /i 7^ in case 2. In the admissible coordinates, up to multiplication of F by —1, case 1 
is equivalent to the condition ac > 0, case 2 is equivalent to the condition ac < and, finally, case 3 is 
equivalent to the condition ac = 0. 

We now consider all three cases. 
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3.2.1 Case 1: ac > 0. 

Without loss of generality we assume a > 0, c > 0. Consider the coordinates (IM)) . In these coordinates 
a = 1, c = 1 and equations (|23[) have the following simple form. 

(fb) y + 2f x = 0, 

(fb) x + 2f y = 0. W 
This system can be solved. Indeed, it is equivalent to 

(fb + 2f) x + (fb + 2f) y = 0, 



(fb-2f) x ~(fb-2f) y = 0, (27) 
which after the (non-admissible) change of coordinates x new ~ x + y, y new = x — y, has the form 

' (fb + 2f) x = 0, 

(fb-2f)y = 0, W 

implying fb+2f = Y(y), fb-2f = X(x). Thus, 

yfa)-^) X(a;)+y(y) 
7 4 ' F(y) 

Finally, in the new coordinates, the metric and the integral have (up to a possible multiplication by a 
constant) the form 

(X-Y)(dx 2 -dy 2 ), (29) 
2 \Px X( x} -Y(y)(Px Py)+Py)- X(x)-Y(v) • ( ' 



3.2.2 Case 2: ac < 0. 



Without loss of generality we can assume a > 0, c < 0. Consider the normal coordinates ([24]) . In these 
coordinates a = 1, c = — 1 and equations (|23|) have the following simple form. 

(fb) y + 2f x = 0, 



(fb) x -2f y = 0. (31) 

We see that these equations are the Cauchy-Riemann conditions for the complex-valued function fb + 2if. 
Thus, for an appropriate holomorphic function h = h(x + iy) we have fb = 8t(h), 2/ = SJ(/i). 

Finally, in a certain coordinate system, the metric and the integral are (up to possible multiplication by 
constants) 

S(h)dxdy and p 2 x - p 2 y + 2^-p xPy (32) 

3.2.3 Case 3: ac = 0. 

Without loss of generality we can assume a > 0, c = 0. Consider admissible coordinates x, y, such that x 
is the normal coordinate from ([24l . In these coordinates a = 1, c = 0, and the equations ((23|) have the 
following simple form. 

' (/&)„ + 2/ x = 0, 



(/&)* - 0. (33) 

This system can be solved. Indeed, the second equation implies fb = —Y(y). Substituting this in the first 
equation we obtain Y 1 = 2f x implying 

f=7zY'(y) + Y(y) and b = . 

2 %Y'{y)+Y{y) 
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Finally, the metric and the integral are 

Y(y) + ^Y'(y))dxdy and p\- Y{y) p x p y . (34) 
2 1 Y(y) + fY'{y) 

Moreover, by the change y new = (3{y id), equations (|34|) will be simply transformed to: 

(?(y)[3' + X -Y' (y)) dxdy and p* x - r^—. . PxPy ■ (35) 

Thus, by putting /3(y) = J^ g Y^dt, we can make the metric and the integral to be 

Y(y) 

PxPy ■ 



(l + ^Y'(y))dxdy and & - ^0 



'{y) 1 

Moreover, after the coordinate change x new = and multiplication of the metric by i, the metric and 
the integral have the form from Theorem [1] 

{l + xY'(y))dxdy and p 2 x - 2 ^^— PxPy . (36) 

Theorem [T] is proved. 

Remark 14. Let us note that if dY ^ 0, then we can take Y as the coordinate y. Then, the metric and the 
integral (f34|) will have the form (see also Remark [4]) 

Y (y) - - ) dxdy and p 2 x + ^— — -p x p y . (37) 



4 Conclusions 

We have discussed integrable geodesic flows of pseudo-Riemannian metrics on 2-dimensional manifolds con- 
structing (Theorem [1]) local normal forms of such metrics. The normal forms are of three types: Liouville 
(the analogous of the Riemannian case) , Complex-Liouville and Jordan-block. We have shown that these 
metrics, in analogy with the Riemannian case, admit geodesically equivalent metrics, can be used to con- 
struct a large family of natural systems admitting integrals quadratic in momenta, that these natural systems 
are integrable by quadratures and that the integrability of such systems can be generalized to the quantum 
setting. A natural further step in this field would be to understand what is the structure of the quadratic 
integral in the case in which the manifold is closed (the Riemannian case is done in [4, 8l [2T1 [22] ). 
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